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The Yale/SMARTS optical-near-IR monitoring program has followed the variations in emission of the Fermi- 
LAT monitored blazars in the southern sky with closely spaced observations since 2008. We report the discovery 
of an optical-near-IR (OIR) outburst with no accompanying 7-rays in the blazar PKS 0208-512, one of the 
targets of this program. While the source undergoes three outbursts of 1 mag or more at OIR wavelengths 
lasting for >3 months during 2008-2011, only interval 1 and 3 have corresponding bright phases in GeV energies 
lasting >1 month. The OIR outburst during interval 2 is comparable in brightness and temporal extent to 
the OIR flares during intervals 1 and 3 which do have 7-ray counterparts. 7-ray and OIR variability are very 
well-correlated in most cases in the Fermi blazars and the lack of correlation in this case is anomalous. By 
analyzing the 7-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the 
location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact 
and the magnetic field strength is higher, and the bulk Lorentz factor of the material in the jet is smaller. These 
result in a much lower Compton dominance and no observable 7-ray outburst during interval 2. 



1. Introduction 



Fermi Gamma-Ray Space Telescope with its clever 
observing strategy of scanning the sky every three 
hours and supporting multi-wavelength monitoring by 
numerous research groups has been ideal to obtain 
high-quality variability information of a large sample 
of blazars with unprecedented richness. The analy- 
sis of the resulting data have have gener ated impor- 
tant r ecent studies of blazar jets (e.g., lAbdo et"al~l 
l2010allbh . The bright source list from the Fermi 2-yr 
catalog (2FGL) contains more than 1000 active galac- 
tic nuclei, most of which are blazars ( Ackermann et all 
l2012t iNolan et al1l2012h . The number of sources is 
large and the variability timescale in many blazars 
is similar at 7-ray and optical-near-IR (OIR) bands. 
Therefore, one should ideally design OIR monitoring 
programs such that the sampling frequency of a source 
at OIR and 7-ray energies are similar while maintain- 
ing a regular cadence for many sources. Given a con- 
stant amount of observing time at the OIR frequencies 
we should arrange the cadence so that sources with 
enough 7-ray flux to be detected with an integration 
time of N days should be observed every N days. This 



means we should increase the cadence of 7-ray-bright 
sources and decrease that of the quiescent ones. 

Here we report one successful example of optical- 
near-IR monitoring with such flexible cadence, 
namely, the Y ale/SMARTS optical-near IR monitor- 



ing p rogram' (|Bonning et al.l I20T2I ; IChatteriee et al.l 
2012). Among other Fermi-LAT monitored blazars 
in the southern sky, we have followed the variations 
in emissi on of PKS 0208-5 12, a blazar at redshift 
z=1.003 (|Healev et al.ll2008h . It was discovered in the 
Parkes Survey of Radio Sour ces dBolton et al.l 1964) 



and was detected regularly ([Thompson et al.l 



1995) 



by EGRET on board the Compton Gamma Ray Ob- 
servatory (CGRO). Most recently, multi- wavelength 
observations o f its la rge-scale jet were presented in 
iPerhnan et ail (|201lh . PKS 0208-512 was originally 
classified as a BL Lac ob ject based on the equiv- 
alent width of Mgll li n e faealev et alJl2008t ). Re- 
cently, iGhisellini et al.l ( 201ll ) pointed out that its 
spectral energy distribution (SED) resembles that of 
FSRQs and its broad lines are very luminous but 
are overwhelmed by the bri ghter continuum luminos - 
ity (L Mg ii ~10 44 erg s" 1 ; [Scarpa fc Falomol Il997l ). 



'Current address: Department of Physics and Astronomy 3905, 
University of Wyoming, 1000 East University, Laramie, WY 
82071, USA; rchatter@uwyo.edu 

t Current address: Quest University Canada, 3200 University 
Boulevard Squamish, BC V8B 0N8, Canada. 



Hence, it should be classified as an FSRQ. 

PKS 0208-512 was detected by the Large Area Tele- 
scope (LAT) on board Fermi at a level higher than 
5xl0 -7 ph cm -2 sec -1 with a time-binning of 1 day 



1 http: / / www.astro.yale.edu/smarts / fermi / 
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Figure 1: Variation of 0.1-300 GeV 7-ray flux from Fermi-LAT in daily and monthly bins, X-ray flux and hardness 
ratio from Swift-XRT, B-, R— and J-band flux densities, and B — J color (from ANDICAM on SMARTS 1.3 m 
telescope) of PKS 0208-512 from 2008 August to 2012 February. The gray shaded regions indicate three intervals 
during which the source underwent long-term and powerful outbursts at optical-near IR wavelengths. While the GeV 
flux increased during intervals 1 and 3, the source was detected at a very low 7-ray flux level at only one monthly bin 
during the second interval. 



regularly from 2008 September to 2008 December and 
from 2011 April to 2011 September. It was too close to 
the Sun at the start of the latter brightening phase. 
We obtained well-sampled optical and near-infrared 
light curves starting 2011 June, when it became ob- 
servable from CTIO, Chile as well as the entirety of 
the interval from 2008 September to 2008 December. 
In addition, we monitored this blazar with one ob- 
servation per night during another OIR bright phase 
from 2009 November to 2010 January during which it 
was quiescent in 7-rays. This OIR outburst is com- 
parable in brightness and temporal extent to the OIR 
flares in the intervals mentioned above during 2008 
and 2011 which do have 7-ray counterparts. In this 
Letter, we analyze the 7-ray and optical light curves 
to investigate, in particular, the physical mechanism 



that can produce this optical-near IR-only outburst. 

In Section 2, we present the observations and data 
reduction procedures. In Section 3, we describe the 
results. We discuss possible scenarios that can explain 
the nature of the observed 7-ray/OIR variability in 
Section 4. 



2. Data 

2.1 . GeV Data 

We derived the 0.1-300 GeV 7-ray flux of PKS 
0208-512 by analyzing data from Fermi-LAT using 
the standard Fermi-LAT Science Tools software pack- 
age (version v9r27pl). We analyzed a Region of In- 
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Figure 2: 0.1-300 GeV 7-ray flux (red open circles) and 
optical R— band flux density (blue filled circles) of PKS 
0208-512 during intervals 1 (top), 2 (middle) and 3 
(bottom). These GeV light curve segments are binned 
over 1-day intervals and include flux values for detections 
with TS > 4.0 (equivalent to ~2cr). The left and right 
hand vertical axes denote the units of the GeV and 
R-band flux, respectively. It is evident that the GeV and 
R— band variability are remarkably well-correlated during 
intervals 1 and 3. 



terest of 15° in radius, centered at the position of 
PKS 0208-512, using the maximum likelihood algo- 
rithm implemented in gtlike, modeling sources with a 
simple power law. We included all sources within 15° 
of PKS 0208-512, extracted from the Fermi 2-yr cat- 
alog (2FGL), with their normalizations kept free and 
spectral indices fixed to their catalog values. We use 
the currently recommended P7SOURCE_V6 set of the 
instrument response functions, Galactic diffuse back- 
ground model, and isotropic background model. In 
the Fermi LAT Second Source Catalog, PKS 0208- 
512 has been modeled with a simple power-law of in- 
dex 2.39±0.04 and its curve significance is 1.6. This 
value of the curve significance is much lower than 16.0 
which indicates switching to a curved spectrum such 
as log parabola will not improve the fit significantly. 
Hence, we fix the spectral index of PKS 0208-512 at 
2.39 in the model hie used in the gtlike analysis. 



The detection significance of PKS 0208-512 in 
2FGL, which integrates the first 24 months of data, 
is 35.8cr. We calculate the fluxes of PKS 0208-512 
from 1-day integrations during 2008 September to 
2012 February, with a detection criterion that the 
maximum-likelihood test statistic (TS) exceed 25.0. 
A TS value of 25.0 is roughly comparable to a 5<r de- 
tection. This light curve is shown in the top panel 
of Figure [TJ The GeV light curve is more sparsely 
sampled than that in the optical-near IR frequencies. 
This is because the blazar was not detectable at TS 
> 25.0 level with 1-day integrations for long periods 
during the 3.5 yr time interval considered here. To in- 
vestigate the nature of GeV variability of this blazar 
during intervals when it was not detected in 1-day 
bins, we also calculate the GeV fluxes from 30-day in- 
tegrations during 2008 September to 2012 February. 
We show the monthly light curve in the second panel 
of Figure [JJ The gaps in the monthly light curve con- 
sist of bins when PKS 0208-512 was not detected at 
TS > 25.0 level even with a 30-day integration. 

We identify three intervals during which PKS 0208- 
512 undergoes outbursts at OIR wavelengths lasting 
for >3 months with date ranges i) MJD 54690-54856 
(2008 August to 2009 January) ii) MJD 55087-55233 
(2009 September to 2010 February) and in) MJD 
55712-55893 (2011 May to 2011 November). These in- 
tervals are shown as the gray shaded regions in Figure 
[TJ These intervals are defined such that i) they con- 
tain a steady rise of flux by 1.3 magnitudes or more ii) 
they contain the corresponding decaying branch down 
to the "quiescent" level at which the rise started hi) 
in the cases where another steady rise of ~ 0.5 magni- 
tude or larger started before the flux level decreased 
to the "quiescent" level, then the interval is cut off at 
the start of the said rise and iv) the length of each 
interval is 2 months or more. While defining "steady" 
rise or decay, we ignore small fluctuations in the ris- 
ing or decaying branch which are less than 0.5 mag- 
nitudes and/or less than a month long because our 
goal in this work is to investigate the intervals con- 
taining longer-term and more powerful outbursts. We 
include two large OIR outbursts together in interval 
3 because they were very close in time and only the 
decaying branch of the first flare was fully sampled 
in OIR frequencies. We assume that for our purpose 
counting these two subsequent outbursts into one in- 
terval simplifies the interval dehnitions without any 
loss of information. We test this in subsequent analy- 
sis. 

In contrast to the three intervals containing large 
OIR outbursts, the source undergoes bright phases in 
GeV energies lasting >1 month only during intervals 1 
and 3. While during intervals 1 and 3, the blazar was 
detected at numerous daily bins and all monthly bins, 
it was detected in no daily bin and only one monthly 
bin during interval 2. The average GeV flux during 
the monthly bin was (5.1 ± 1.4) xlO -8 ph cm -2 sec -1 . 
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Figure 3: Discrete cross-correlation function (DCCF) of the 7-ray and _R-band light curves of PKS 0208-512 during 
intervals 1 and 3. The time delay is defined as positive if the R— band variations lead those at 7-ray energies. During 
both intervals, the 7-ray and R— band variations are very strongly correlated with zero time lag. 



This is more than an order of magnitude fainter than 
that in intervals 1 and 3. To better compare optical 
and 7-ray variations, we also obtained a GeV light 
curve of this blazar with a detection criterion of TS 
> 4.0 (comparable to a 2a detection). We plot this 
light curve for all three intervals along with R— band 
variability in Figure [2] 

2.2. Optical-Near-IR and X-Ray Data 

All of the measurements in B-, R-, and J— band 
are from the ANDICAM instrument on SMARTS 
1.3m telescope located at CTIO, Chile. ANDICAM 
is a dual-channel imager with a dichroic that feeds an 
optical CCD and an IR imager, which can obtain si- 
multaneous data at one optical and one near-IR band. 
For details of data acqui sition, calibration an d data 
reduction procedures, see iBonning et all (|2012h . 

We obtained the X-ray flux and hardness ratio dur- 
ing the same interval from the Swift-XRT monitoring 
program of Fermi-L AT sources of interest 2 . Most of 
the observations were carried out with exposure times 
of ~1 ks. The mean count rate was 0.07 cts/s in 
0.3-10.0 keV energy range and the mean background 
count rate was 0.0007 cts/s during this entire period. 
As in the optical, X-ray variability was present during 
all three interva ls. We use the Sw ift-XKT data prod- 
ucts generator 3 ( Evans et al.ll2009h to derive the aver- 
age photon index during the three relevant intervals. 
The photon index stays between 1.58 and 1.75, indi- 
cating no significant difference between intervals. The 
gaps in the X-ray light curve denotes intervals during 
which PKS 0208-512 was not observed by Swift. B-, 
R-, and J— band light curves, and the variation of the 
B — J color, X-ray flux and hardness ratio are shown 
in Figure [TJ 



2.3. Parsec-Scale Structure of the Jet 
from VLBI data 

Very Long Baseline Interferometry (VLBI) imaging 
is the only direct way to study the dynamics of blazar 
jets at parsec scale. Contemporaneous VLBI moni- 
toring along with photometric monitoring at multiple 
wave bands can provide important clues about the lo- 
cation and mechanism of relevant emission processes 
in blazar jets. Time variability studies in the EGRET 
era involving GeV light curves and Very Long Base- 
line Array (VLBA) monitoring have indicated a con- 
nection between the "ejection" of new radio knots in 
the pc- scale jet of blazars and outbursts at 7-ray en- 
ergie s ( Jorstad et al.l l200ll lLahteenmaki fc Valtaoial 
2003). Ejection is defined as the separation of a radio 
knot from the "core". The core may be a standing 
shock seen as a (presumed stationary) bright spot in 
the pc-s cale jets of bl a zars i n VLBA images. More 
recen tly, lAgudo et al.1 (|201ll ): iMarscher et al ] (|2008l . 
2010) have inferred that the 7-rays in some blazars are 
produced downstream of the VLBI core (which lies >1 
pc from the central super-massive black hole) analyz- 
ing a combination of time-dependent multi- waveband 
flux and linear polarization observations, and sub- 
milliarcsecond-scale polarimetric VLBI images at 7 
mm. 

The blazar PKS 0208-512 is a source in the southern 
sky and it has been observed roughly once every six 
months since the end of 2007 as part of the Tracking 
Active Galactic Nuclei with Austral Milliarsecond In- 



2 http: / /www. swift. psu.edu/monitoring/ 
3 http: / /www. swift. ac.uk/user_objects/docs.php 



terfer ometry (TANAMI) campaign (jBlanchard et al.l 
2012). The VLBI observations are typically made 
with all Australian LBA antennas (the Australia Tele- 
scope Compact Array or ATCA, Parkes, Mopra, Ce- 
duna and Hobart) and often other southern hemi- 
sphere antennas including the Tidbinbilla 70m, TIGO, 
O'Higgins and Hartebeesthoek. The images from this 
monitoring program shows that a new jet component 
may have been ejected ^6—10 months before both 
large 7-ray outbursts (intervals 1 and 3) but not be- 
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fore or during interval 2. Assuming an uncertainty of 
~3 months because of the sampling rate of once every 
6 months in this program, it may imply a connection 
between an increase in GeV emission and ejection of 
new components in the pc-scale jet similar to that seen 
in other blazars, e.g., PKS 1510-089, OJ 287, and BL 
Lac. 



3. Results 

3.1. 7-Ray/Optical Correlation in 
Intervals 1 and 3 

From Figure [2j it is evident that the GeV and 
R— band variability are remarkably well-correlated 
during intervals 1 and 3. To confirm this, we cross- 
correlate the 7-ray and R— band light curves us- 
ing the discrete cross- correlation function (DCCF; 
lEdelson fe Kroliklll988l) . A similar cross-correlation 
analysis for interval 2 was not possible due to infre- 
quent detection even at TS > 4.0 level. As shown 
in Figure |3l the 7-ray and i?-band variability in this 
blazar during both intervals are very strongly corre- 
lated with zero time lag. This is consistent with simi- 
lar correlated 7-ray and optical/near-in frared variabil- 
ity s e en in many other blazars (e.g., iBonning et al.1 
l2009l . 120121) . The degree of correlation and lack of 
a significant lag can be explained by the standard 
leptonic scenario in which the 7-rays and OIR emis- 
sion are generated by the same relativistic electrons 
in the jet through inverse- Compton and synchrotron 
processes, respectively. The source of the seed pho- 
tons that are being scattered may be the synchrotron 
photons generated within the jet, in which case it 
is termed synchrotron se l f-Compton (SSC) pro cess 
(jChiang fc Bottcherl [200l iMaraschi et al.1 11992T ). or 
from outside the jet (radiation from broad emission 
line region, accretion disk, or dus ty torus) , termed ex- 
terna l Compton or EC p r ocess (ICoppi fc Aharonianl 
I1999L iDermer et all 120091 : ISikora et all Il994h . Both 
SSC and EC scenarios predict strong correlation with 
zero time lag. 



3.2. Absence of Correlation in Interval 2 

As shown in Figure [T] the source underwent an 
optical-near IR outburst during interval 2. The in- 
crease of brightness was comparable to that in inter- 
val 1 and the B — J color variations (bottom panel of 
Figure [T]) were similar to those in intervals 1 and 3. 
But unlike those flares, in interval 2, there was no sig- 
nificant increase in the GeV flux and the average GeV 
flux was smaller than average flux values of intervals 
1 and 3 by at least an order of magnitude. 

Note that there is no significant difference in the na- 
ture of the GeV/R— band correlation in the two con- 



secutive flares during interval 3 (Figure [2]). This is 
consistent with including them together as part of the 
same interval. 



4. Discussion 

There is one scenario in which a change at optical- 
near-IR wavelengths is not accompanied by a change 
in the GeV energies, namely, the optical synchrotron 
emission is due to a change in the magnetic field 
only, i.e., its magnitude and/or direction, while the 
GeV emission is due to a temporarily steady exter- 
nal Compton process. Inverse-Compton emission de- 
pends on the total number of emitting electrons (N e ), 
Doppler factor (<5), and the number of seed photons 
available for scattering while synchrotron emission de- 
pends on N , S, as well as the magnetic field and the 
viewing angle in the form Bsin0, where viewing angle 
(</)) is defined as the angle between the direction of the 
magnetic field vector and the line of sight. Variability 
in the observed 7-ray and OIR emission can be gener- 
ated by changes in any one or a combination of these 
parameters. Hence, we conclude that the OIR out- 
burst during interval 2 might be caused by a change in 
the magnetic field in the emitting region without any 
change in the other two parameters described above. 
The X-ray emission from F SRQs is often 



domi n ated by the SSC p r ocess dChatteriee et all 
2008 ; iHartrnan et all l200lt iMukheriee et all Il999f 
Sikora et al.ll200ll) which means that changes in mag- 
netic field should lead to X-ray variability at a level 
similar to those changes. This is roughly consistent 
with the relative variability in optical and X-ray 
flux seen here. We note that the optical synchrotron 
emission and the SSC X-rays are produced by differ- 
ent parts of the electron distribution. Assuming a 
B field of 1 Gauss, synchrotron emission at optical 
wavelengths is generated by electrons of Lorentz 
factor 7 ~ 10 4 while SSC X-rays are produced by 7 
~ 10 2 assuming average seed photons to be at the IR 
wavelength range where the synchrotron peak occurs 
for most FSRQs. The difference in the OIR and 
X-ray variability, if any, in 0208-512 during interval 2 
can be attributed to this. 

In another possible scenario, the large difference 
in the GeV/OIR ratio between intervals 1 and 3 
vs interval 2 may be caused by the location of 
the outbursts. The acceleration and collimation of 
magnetohydrodynami cally launched jets occur over 
an extended region dMcKinnev fc Blandfordl [2 009; 



Meier. Koide. feUchidal 120001 : IVlahakis fc Konigll 
2004), possibly as high as ~10 4 times the 
Schwarszchild radius. Hence, an outburst tak- 
ing place closer to the black hole may be associated 
with a smaller bulk Lorentz factor (r). In contrast the 
corresponding magnetic field and particle density will 
be larger due to its compactness. It has been shown 



eConf C121028 



4 th Fermi Symposium : Monterey, CA : 28 Oct- 2 Nov 2012 



that the EC/synchrotron ratio can be dramatically 
different depending on the location where a given 
outburst takes place, with EC becoming relatively 
more do minant for outbursts occurring farther down 
the jet ( Katarzvhski fc Ghisellinill2007l ). The appar- 
ent seed photon field seen by the emitting electrons 
is boosted by a factor of T 2 causing a significant 
dependence of the EC /synchrotron ratio on T. If 
the outburst during interval 2 is generated closer 
to the BH where Y is smaller while the outbursts 
in intervals 1 and 3 are produced farther down the 
jet where Y is larger, the optical-near IR emission 
will be significantly more dominating than that in 
the GeV bands during interval 2 than the other 
two. A factor of two decrease in Y coupled with a 
factor of two increase in B can produce the change of 
EC/synchrotron ratio we see from interval 1 to 2. 

Modeling of truly simultaneous time variable spec- 
tral energy distribution (SED) from radio to 7-rays 
is the holy grail of blazar physics. Such analysis 
for a large sample of blazars is now possible with 
the data available from Fermi and supporting multi- 
wavelength programs. However, detailed spectral 
studies are quite ambitious in their need of spectral 
coverage and computation time. Given a stipulated 
supply of computation time and multi- waveband data, 
suitable events need to be sorted out for detailed 
study. This study is one of the very few cases where 
one of two specific physical scenarios can be identi- 
fied for certain observed events in a blazar jet. These 
events are ideal candidates for detailed modeling of 
the time variable SED and inspired multiple model- 
ing studies which are currently underway (Chen et al., 
Chatterjee et al., both in prep). These may provide 
stronger constrains on the relevant physical parame- 
ters in blazar jets which may lead to a better under- 
standing of the physical parameters and processes in 
the jet very close to the black hole which is one of the 
fundamental goals of studying blazars. 
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